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Abstract 

Perturbative contributions to single-beam two-photon transition rates may be divided into two 
types. The first, involving low-energy intermediate states, require a high-order perturbation treat- 
ment, or an exact diagonalization. The other, involving high energy intermediate states, only 
require a low-order perturbation treatment. We show how to partition the effective transition 
operator into two terms, corresponding to these two types, in such a way that a many-body per- 
turbation expansion may be generated that obeys the linked cluster theorem and has a simple 
diagrammatic representation. 
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I. INTRODUCTION 



Two-photon laser spectroscopy is an important complementary technique to linear spec- 
troscopy because it has a different parity selection rule, allows access to higher energy states, 
and has a greater variety of possible polarization choices than linear spectroscopy-* The cal- 
culation of two-photon absorption for strongly correlated many-electron systems, especially 
for rare earth ions in solids, dated back to Axe's work in 1964.— However, few quantita- 
tive measurements of two-photon absorption were made until the early 1980's.— Extensive 
measurements have been carried out thereafter . 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ' 11 Most of these measurements 
cannot be explained by Axe's lowest-order calculations, and as a result, many calculations, 
some using perturbation theor y 5 i 12 i 13 i 14 i 15 i 16 i 17 and some using full calculations in truncated 
spaces , 18 ! 19 have been carried out to interpret the experimental results. Recently, we have 
been able to explain the puzzling two-photon absorption intensities and polarization de- 
pendences of Gd 3+ :LaF 3 and Eu 2+ :CaF 2 by full calculations in a truncated 4f N + 4f N ~ 1 5d 
space . 19 ! 20 However, these calculations cannot explain the two-photon absorption intensi- 
ties of Sm 3+ , Eu 3+ and Tb 3+ doped in elpasolites . 21 ! 22 ! 23 It appears that contributions from 
high-order perturbations or high energy intermediate states must be considered.— 

A systematic way to calculate properties of many-body systems is using effective 
Hamiltonians and operators. These techniques have been extensively developed in the 
literature . 25 ! 26 ! 27 ! 28 ! 29 ! 30 The basic idea is to transform the exact time-independent (usually 
many-body) Hamiltonian H into an effective Hamiltonian H e g acting on a restricted model 
space of manageable dimension. The exact eigenvalues and model space eigenvectors (not 
the exact eigenvectors) can be obtained by diagonalizing H e g. For a time-independent oper- 
ator O, such as a transition moment operator, an effective operator O e s may be introduced 
that gives the same matrix elements between the model space eigenvectors of H e g as those 
of the original operator O between the corresponding true eigenvectors of if.— Although 
the forms of if e ff and O e s are generally more complicated than, respectively, H and O, the 
calculations based on H e ff and O e s have many advantages over variational and other direct 
calculations based on H and O, such as smaller bases, less calculation effort, order by order 
approximations, and the calculation of all eigenvalues and transition matrix elements simul- 
taneously for a multi-dimensional model space. More details on effective operators can be 
found in a recent paper of Duan and Reid^ and references therein. 
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Many-body perturbation theory (MBPT) may be used to expand effective Hamiltoni- 
ans and operators order by order. The calculations are often represented by Goldstone 
diagrams . 31 ! 32 The linked cluster theore m 25 i 32 i 33 implies that disconnected diagrams cancel 
for effective Hamiltonians and effective operators, provided that the model space has been 
carefully chosen and the model states have been properly orthogonalized. This cancellation 
reduces the number of high order diagrams greatly and ensures size consistency. The linked 
cluster theorem also holds for one-photon transition operators . 25 ' 29 However, the application 
of MBPT to two-photon transitions is more difficult than the case of one-photon transitions, 
in that there are energy denominators in two-photon transition operators that contain both 
the photon energies and exact electronic energies.— >2i Also, the intermediate states can be 
any eigenstates of the system, including states in the model spaces containing initial and 
final states and other low excited states, making the energy denominators change drastically 
and makes perturbative expansions for these intermediate states impossible. 

In this paper, we explore the effective operator method for two-photon transition calcu- 
lations by combining exact calculations in a truncated spac o 1 ^ 19 with perturbative methods 
for the rest of the states . 171 In section II we review the basic formalism for effective operator 
methods; In section III the partition of two-photon transition operator is given in detail; 
Section IV presents the perturbation expansion that may be suitable for diagram repre- 
sentation and applying linked cluster theorem. The diagram representation and diagram 
evaluation rules and linked cluster theorem themselves for effective two-photon transition 
operator are highly nontrivial and will be presented in a followed paper. 

II. BASIC FORMALISM 

Most of the formalism required has been treated in detail in a monograph by Lindgren 
and Morrison 32 and more recently summarized by Killingbeck and Jolicard.^ The concept of 
biorthogonal model space eigenvectors has been summarized by Duan and Reid.- 2 ^ Here shall 
only give a brief description of the formalism necessarily in the presentation that follows. 

The time-independent Hamiltonian H is written as the sum of a model Hamiltonian H 
and a perturbation V 

H = H + V. (1) 
Usually H is chosen in such a way that its eigenvalues and eigenvectors can be obtained 
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more easily than for H. For example, when H is the Hamiltonian for a many-body system, 
Hq is usually chosen such that each particle moves independently in the average field of 
other particles and V denotes the remainder of H. 

A complete set of orthonormal eigenvectors {|a)} and corresponding eigenvalues {Eq} 
for H are assumed to be available 

H \a) = Eg\a), (2) 
(a\(3) = V (3) 

A general model space Po, often referred to as a quasi-degenerate model space, is defined as 
the space spanned by d successive eigenvectors of H Q (not necessarily strictly degenerate). 
The remaining part of the Hilbert space is called the orthogonal space Qq. Here we denote 
the associated projection operators also as Po and Qq: 

d 

P = X>)H (4) 

Of=l 

Qo = J2\a)(a\ = l-P (5) 

a>d 

It has been shown that d eigenvectors of the full Hamiltonian \^ a ) (a = 1, ■ • ■ , d) can usually 
be projected into the model space as d linearly independent functions |\I/q) in P . 27 ' 35 The 
wave operator Qp is defined as an operator that transforms all the d model functions back 
into the corresponding exact eigenvectors. 

|*g [ > = Po|* a > (6) 
\^ a ) = Q P \%). (7) 

Note that |^o) s are n °t necessarily eigenstates of H , but a linear combination of those 
eigenstates of H in P . We shall call the space spanned by the d exact eigenvectors P. The 
wave operator Qp satisfied the ordinary Bloch equation^ 

[n Pj h q ] = (vn P - n P p vn P )p . (8) 

Usually the |\l/g) are not orthogonal but are chosen to be normalized to unity. As a 
consequence is not normalized to unity. Denote 

= N% p . (9) 
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An effective Hamiltonian H e s acting on the model space, which gives the d exact eigen- 
values and model eigenvectors upon diagonalizing, can now be defined. Its form and eigen- 
equation are 

H eS = P Hn P P = P H P + PoVtlpPo, (10) 
H cS \^) k = E a \%) k . (11) 

Instead of calculating Qp directly from (J8|i. the effective Hamiltonian are usually calculated 
via perturbation theory or phenomenological method and then diagonalized to give eigen- 
values E a and eigenvectors 1^0%, It is straightforward to find from a set of vectors b(^o \ in 
the model space such that 

b (V%\*g) k = <W (12) 
|*J>* 6<*?| = Po (13) 

It is straightforward to show that 

bi^m+np)- 1 e P , (14) 
b (^|(fi+fip)- 1 fi+fi P |^) fc = 5 aP , (15) 
„(%\H eS = E a b (%\, (16) 

which together with Eq. Eland Eq. El show that 

b (^\(Q + pQp)- 1 Q + P = N- 2 (^\, (17) 

p = J2 npPointnp)- 1 ^. (is) 

a 

The transition matrix element of an operator O between states (a G A) and (\E^| 
((3 e B) is 

° aP ~ N a N/3 ■ (19) 



Defining 



it can be shown that 



O cfr = {Qp^ttp^Qp-OQp, (20) 



h (*%\0«\*g) k = ^0 Q/3 , (21) 
O a/3 = ( b (^\O eS \<)k ,(^10^)1)-^. (22) 
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III. PARTITION OF EFFECTIVE TRANSITION OPERATORS FOR SINGLE 
BEAM TWO-PHOTON ABSORPTION 



The perturbative expansion of a general two-beam two-photon transition rate has been 
discussed by Duan and Reid^. Here we develop the perturbative expansion specialized for 
the single-beam case. We will see that the calculation can be greatly simplified. 

The line strength is proportional to the square modulus of following expression 

^ _ (f\Q\k)(k\Q\i) 
Ei + lo — Ek 

where |/) (/ G F), \i) (i G /) and \k)(k G F U / U K) are normalized exact eigenstates of 
the systems, Ef, E iy E k are the corresponding exact eigenvalues, F, I are the set of final 
and initial states of the transition. K is the set of states of the system not included in 
/ and F . The intermediate eigenstates k can be any eigenstates of the system, including 
eigenstates in the sets of initial and final states. The calculation of T/j can be divided into 
two terms, a term 7\ with "small" denominators, where k G J U F, and a term T 2 with 
"large" denominators, where k G K. The operator can be formally written as 

T = T 1 + T 2 (25) 

t 1 = El/M E , E ° f E k °;: E (26) 

f,i keiuF^Ef + Ei)/2- F k 

- ( m + w)i2-H < 27 > 

T 2 = El/X'l E (E E (28) 

= v+A/2- g °- (29) 

Note that H F and H 1 are actually acting on the transition final (on the leftmost) and 
initial states (on the rightmost) and the equalities in (}27j) and (|29|) above are for notational 
convenience. 

The effective operator for T 1 is 

f,i m&FUl [Ef + Fi)/2- F m 
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Note that in this expression we use the exact eigenvalues and model space eigenvectors. The 
matrix elements of the effective operators are between states in spaces Ppo and Pjq (which 
may be the same). The matrix elements may be calculated from the effective Hamiltonian 
H e g and effective operator O e g. 

The effective operator for T 2 can be formally written as 

T c 2 ff = (p$n F )-Hi$on K p% F J ^ ^(n+n^-^+onz. (31) 

v cff ' cff / / ^ei 

Once again, H^ s and H^ act on the bra and ket model space respectively. 

The space K includes all other states of the system than those limited number of states 
in / and F. It is usually of infinite dimension and the calculation of H^, O e g(F, K) and 
O e s(K, I) is usually impractical or at least very tedious. Perturbative expansions giving in 
the following section can be used to calculate T e 2 ff by an order-by-order approximation. 



IV. PERTURBATIVE EXPANSION 



Perturbative expansions of H e Q and O e g have been discussed in, for example, Refs. |32 



and |29| by applying the Bloch equation iteratively. Here we expand T e 2 ff by perturbation 
theory to avoid direct calculation of H cS and O cS in model space K, which is usually of 
infinite dimension. 
Defining 

S = l/[(H* + Hl)/2-H*), (32) 

= Kff-^Kff + Kff), (33) 

where 

V& = PtfWVfc (34) 

Hi = PzoHoPzo (35) 

act on the transition initial, final and intermediate states for £ = I, F, K respectively. The 
energy denominator can be expanded as follows 

mTHSsj^HS =s ^ VS) ' (36) 

= S + SVj'S - iv £S ! _ I S VJ, + ■ • • (37) 



Using the Bloch equation, VLk can be expanded as follows 



n K = p Ko + R K (vn K - n K vn K ) (38) 

= P Ko + RkVP K0 + R K VR K VP K0 - R 2 VP' K0 VP K0 , (39) 



where 



Rk ~ H«-H^- (40) 
Using the above expressions, the matrix elements of the zeroth and first-order of T 2 S between 
eigenstates of H , \ f) (/ G Pfo) an d |z) {i G P/o) are: 



fce p K0 (E f0 + E t0 )/2- E k0 



(f\T* s>1 \i) 

+ E 




„ (/|y|Q(i|0|fc)(fc|o|t) 

zgq fo ~~ ^o)[(P/o + Pio)/2 - Pfco] 
(f\0\k)(k\Q\l)(l\V\z) 



+ E E 
+ E E 



(/|0|/)</|V|*)<*|0|i) + (/|0|*><*|V|J>(J|0|i) 



(Eko — Ei )[(Efo + E i0 )/2 — E k0 ] (E k0 — Ei )[(Ef + E i0 )/2 — E k0 ] 

(fioihjihivi^ihioii) 



kl& P K0 k2ePl(0 [(Efo + E i0 )/2 - E kl a\[(Ef + E i0 )/2 - E k20 ] 

_1 j2 V (/|V|f)(/lQ|fc)(fc|OK) 
2 fe p F0 fee p K0 [{Ef + E i0 )/2 - S fc0 ][(-B/'o + E i0 )/2 - E k 

_ i ^ ^ (/|Q|fc)(fc|Q|^(»iy|») , (42) 

2 i, gPffl fc6 p A . [(P/o + £io)/2 - P fco ][(P/o + E vo )/2 - E kQ ] 

where all eigenvectors and energies are for model Hamiltonian Hq and matrix elements are 
between eigenvectors of Hq. The transition rates can then be calculated straightforwardly 
from Eq|221 since the model space eigenvectors are assumed to have already been calculated 
from EqJIll 

Terms of second order or higher in V can also be obtained straightforwardly. There are 
about 20 second order terms but around 100 third order terms. Fortunately, with a suitable 
partition of T into T 1 and T 2 , usually only the zeroth and first order terms of T 2 need to 
be calculated, except when zeroth and first -order terms become zero due to selection rules. 
In such cases the number of nonzero second order terms is often greatly reduced. 
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V. CONCLUSION 



A method to calculate single-beam two-photon absorption transition rates for many- 
electron systems has been developed using effective operator methods together with many- 
body perturbation theory. In this method the contributions to two-photon transition oper- 
ator are partitioned into two terms, one with small drastically varying denominators, which 
is treated by doing an exact calculation in truncated spaces and the other with numerous 
intermediate states and large energy denominators, which is treated systematically with 
many-body perturbation theory Compared to previous methods, the method presented 
here has the accuracy of full calculation for contributions due to drastic-varying low-energy 
intermediate states and the simplicity of low-order many-body perturbation theory for con- 
tributions due to high energy intermediate states. It is also expected that there are linked 
diagram representations for the order-by-order expansion. 
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